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By examining the bound-state spectra of transverse fluc- 
tuations about one-dimensional, spatially localized dark and 
bright soliton wavetrains of the Gross-Pitaevskii equation, it 
is established that the low-temperature ground states of repul- 
sive and attractive quasi-one-dimensional Bose-Einstein con- 
densates are degenerate. In the one-soliton limit, both ground 
states are shown to possess two distinct transverse fluctua- 
tion modes which can couple to the spatial soliton: the first 
involves zero-energy exchange, costing only the soliton shape 
dressing via uniform translation of its centre of mass. The sec- 
ond mode contributes in a negative energy for the repulsive 
case, but positive energy in the attractive case. This unsta- 
bility of the repulsive spatial soliton against localized trans- 
verse fluctuation modes invalidates the Gross-Pitaevskii equa- 
tion for stationary solitonic processes in repulsive ID Bose- 
Einstein-Condensate systems. 



It was known for long time that dilute systems of Bose 
gases can become unstable against various kinds of inter- 
actions. However, the recent observations |jl|-^ of collec- 
tive excitations in Bose-Einstein Condensed(BEC) gases 
motivated pressing renewal of interest to the problem in 
several, yet envisaged theoretical aspects. Solitonic be- 
haviours Q are thus currently being the most wanted 
perspectives in the light of recent experiments 
The framework is the Mean-Field(MF) theory which 
provides acceptable mathematical tools to deal with the 
many-body character of the interacting system. The 
underlying concept sits in the possibility to cast a col- 
lection of microscopic processes to a macroscopic scale. 
At this scale, the order parameter relies on the macro- 
scopic density of the excited atoms then suggesting their 
"quasi-particle"-like properties as they interact with the 
MF(internal) and trap(external) environments. Depend- 
ing on the atomic species, the s-wave scattering length 
for atomic colhsions will be positive [|j or negative 
such that the MP environment is either repulsive (self- 
defocusing effect) or attractive(self-focusing effect) vis- 
a-vis the Bose quasi-particle. The first effect promotes 
collective excitations in terms of "quasi-particle holes" 
with shape preserving macroscopic density. As for the 
second, it gives rise to quasi-particles in form of well local- 
ized hump whose macroscopic density also is shape pre- 
serving. These two different quasi-particle architectures 
are reproductible in the order parameter via the Gross- 



Pitaevskii(GP) equation governing the dynamics of 
the ground-state wavefunction of a one-component BEG 
system. Concerning the true nature of this ground state, 
several questions are still seeking answers but however, 
there are serious indications of a dominant regime where 
the condensate tends to be frozen and thus oscillate 
rather slowly. This suggests the emergence of spatially 
localized collective excitations in the condensate , the 
dynamics of which are well approximated in the MP de- 
scription by a Hatree-Pock treatment pT| ] and which, in 
the contexts of one-dimensional(lD) and quasi-lD(QlD) 
BEG systems, can readily be looked on as stationary 
states p^-[l^ of the GP equation. Regardless of dimen- 
sionality, the GP equation is comparable to a ND non- 
linear Schrodinger(NDNS) equation with the trap-field 
potential playing role of perturbation. The NDNS equa- 
tion(actually the ID version) is popular in the literature 
of optical communication systems jl^ for the richness of 
its soliton solutions: dark and bright solitons of this equa- 
tion are representives of self-defocused and self-focused 
excitation modes, respectively. 

To a viewpoint of the general approach to perturbed 
IDNS equations their equivalence with the GP equa- 
tion means that at least at low enough temperatures, 
the trap does not interfer in the initial stage of the soli- 
ton creation but instead, acts as a drive tightly confin- 
ing the quasi-particle motion in a given direction along 
wich the soliton gets accelerated, increasing its lifetime. 
To this point, it may be instructive to draw attention 
on the fact that most of the currently reported typi- 
cal quasi-particle lifetimes are short enough compared 
to the typical experimental time scales. Therefore to al- 
low direct observations of their oscillations over control- 
lable time periods, it is useful to increase their lifetimes 
by microtrap techniques [0, hollow-blue-detuned laser 
beam traps and so on. On the other side, still in 
the spirit of the perturbation theory fo IDNS equations, 
one should admit that effects of the trap can be advan- 
tageous in stabilizing the soliton wave-shape profile as 
shown in recent attempts assuming nonlinear periodic 
potential traps . Moreover, they can cause soliton 
decays into distinct structures Q] . However these effects 
fall out of our context and will not be considered. Instead 
we assume an homogeneous potential field that allows us 
shifting the zero-point energy up to the desired charac- 
teristic energy scale. 

The recent works on the solitonic properties of BEG sys- 
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terns focused mainly on the pure ID regime(assuming 
that the characteristic length scale of the condensate in 
a given direction is far greater compared to the two oth- 
ers). A priori, a pure ID description is unadequate and 
gives a rather oversimplified view of the physics. In- 
stead the idea of a QID process fits better with exper- 
iments ]T9|-|2^ . One of the relevant peculiarities of the 
2D context is the possible coupling of transverse fluctua- 
tions to the confined ID collective low-frequency oscilla- 
tions in the ground-state phase of the condensate ||2^ . In 
view of the important role played by fluctuations in the 
occurence of a true condensate phase transition, it is sur- 
prising that only little attention has so far been devoded 
to the question. The few previous descriptions of trans- 
verse fluctuations in QID and Q2D BEG systems ||2^ ] 
postulated the interacting gas to be effectively kinemat- 
ically 2D, but undergoing zero-point oscillations in the 
transverse(frozen) direction while the collective motion 
of the condensate is along the confinement. The asserted 
goal of this assumption is quite appealing since it ques- 
tions the possibihty for the condensed ground state to 
be degenerate. However, a preriquisite to this question 
is the ones about the structural (wave shapes) and stabil- 
ity properties of the modes contributing to this degener- 
acy. Current stability theory for ID soliton systems can 
provide relevant insights to these last questions. In this 
respect, we solved the QID GP equation including the 
coupling between the confined ID collective condensate 
excitations and the transverse fiuctuation modes taking 
burn in reponse to the creation of these ID collective ex- 
citations. We found a rather rich bound-state spectrum 
for both repulsive and attractive systems which indicate 
ground states that are indeed degenerate. However, in 
the same time we observed that the current model could 
not be valid for at least one of the two regimes and hence 
need be reformulated. 

In the first step we give a brief, but self-consistent sum- 
mary of the general approach to the soliton stability in 
the presence of phonon fluctuations. This theory is spe- 
cific to the current problem since we are interested in 
periodic spatial solitons 0,^. In the next step, the the- 
ory will be applied to the repulsive and attractive BEG 
systems within the GP equation. Extension to the single- 
soliton limits will reveal itself more interesting since cor- 
respond to most of the current considerations. 
Starting, consider an eigenvalue problem describing the 
scattering of a quasi-particle of finite energy and mass, 
in a ID nonlinear periodic potential: 
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P{k, e) depends on the eigenvalue e associate to the eigen- 
function /(y, k), and sn is the Jacobi Elliptic function of 
real modulus k between and 1. The limits k = 1 and 
n = 1,2 are well known as " Goldstone- Jackiw" eigen- 
value problems HJ]. For K ^ 1 and integer n, equa- 
tion (|l|) turns to a n*''-order Lame equation ||2^] whose 



bound state spectrum has already been discussed [g6 27 
for n = 1 and n = 2. We improve these previous meth- 
ods by introducing more simplified considerations lead- 
ing straightly to the eigenvalues and eigenfunctions. In 
connection with the physical problem to be considered 
below, we will focus on the case n = 2. In general, how- 
ever, for even values of the integer n the bound state 
spectrum will possess 2n + 1 eigenstates corresponding 
to the 2n + 1 possible values of P{k, e). Their wavefunc- 
tions are all real in y and obey the relation: 

/(y, k) = W{y, K)sn'{y)cn''{y)dn'^{y), (s, c, d) = 0, 1. 

(2) 

W{y, k) is a polynomial in sn'^{y) of degree (n — s — c — 
d)/2. Any two combinations of the nonzero values of s, c 
and d will give an eigenfunction of the eigenvalue equa- 
tion provided the total degree of /(y, k) remains equal to 
71. The eigenfunctions for s = c = d = are obtained by 
writing down the polynomial Wijj, k) as: 
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The sum in ^ runs over N{= n/2) such that coeffi- 
cients \2r having indices in between and N are all 
nonzero, and those of indices greather than N are all 
zero. The compatibility relations obtained by replacing 
all the wavefunctions in equation (|l|) form a {2n+lY ma- 
trix equation, the 2n+ 1 eigenvalues of the matrix being 
the 2n + 1 allowed values of P{K,e). Thus constructed, 
the discrete spectrum of the eigenvalue problem (Q) forms 
an orthonormal basis on the associate wavefunctions. 
Let us now see how this theory can help us solving the 
problem of bound fluctuation modes in QID BEG sys- 
tems. Assuming spatially homogeneous potential trap 
/io, the GP equation turns to the perturbed 2DNS equa- 
tion i.e, 
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with ip{f, t) the macroscopic order parameter of the con- 
densate and A the 2D Laplacian. The paramter g is 
the strength of the MP interaction. It is customary to 
write g = ^ where £ is the s-wave scattering length for 
inter atomic collisions and m the mass of atoms. Posi- 
tive and negative values of £ correspond respectively to 
a self-defocusing effect and a self-focusing effect of the 
medium on the BEG quasi-particle. We are interested 
in ID collective excitations undergoing small-amplitude 
shape dressings induced by transverse fuctuations that 
couple to their spatial equilibrium positions ^8 29 1 along 
the confinement. In the coherent regime, the macroscopic 
wavefunction describing such dressed solutions can be 
represented by: 
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Po is the uniform atomic density in the condensed phase, 
is the component of the macroscopic order parame- 
ter describing pure ID coherent quasi-particle excitations 
with ernergy /i and we assume the condensate phase at 
equiUbruim Oq. The function (p'' represents the spatial 
dressing experienced by the wavefunction when trans- 
verse fluctuations p±_ of the amplitude of the order para- 
mater are taken into consideration. The form (|^) of the 
wavefunction solution is appropriate to describe ID con- 
finements in the QID system. Inserting (^) in equation 
and keeping only linear terms in the fluctuations, we 
arrive at the following system of coupled equations: 



2m dx^ 
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Where we define the shifted zero-point energy as Ap — 
p — Po- In the presence of the MF interaction, The order 
of magnitude of this zero-point energy is determined by 
the requirement that in the limit t oo the solution of 
the GP equation (||) effectively describes the true equilib- 
rium state of an interacting Bose gas. That is, neglecting 
all fluctuation terms the wavefunction must be repre- 
sentative of a complex order parameter with amplitude 
and phase which are macroscopically well defined at the 
equilibrium. Therefore it is only with time evolution and 
spatial modulations that the amplitude and phase of the 
order parameter will start fluctuating. In general at low 
dimensions these fluctuations lead to short-range corre- 
lations indicating the formations of topological defects. 
However, by the hierarchy of times scales and character- 
istic length scales involved with fluctuation phenomena 
in BEG systems, these topological defects are rather var- 
ious and relate to both amplitude and phase fluctuations. 
Nevertheless, the relaxation of phase fluctuations gener- 
ally requires longer time such that if one was interested 
in ID short-range processes governed by an order pa- 
rameter varying in time and space like (o), according to 
equation the MF value gpo will constitue an appropri- 
ate energy scale. It turns out that the zero-point energy 
Ap is about this scale. With these considerations, we 
can readily define the characteristic lengths for the 
ID coherent excitations of the order-paramater ampli- 
tude within the GP equation (||) either in terms of the 
zero-point energy, or in terms of MF energy scale. We 
consider the form: 
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With help of the incertainty principle, we can also intro- 
duce a characteristic time scale: 
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As the energies of all the macroscopic processes governed 
by equation are now contrained to the MF scale, it is 
natural to expect topological defects. With respect to the 
two signs of the MF strength g, two possible nonlinear 
topological solutions of this equation are: 
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Parameters k are such that < k < 1 for ( [ll] ) and 
- j^ < K < 1 for (p^. Ap and g are both positive in 

(y) but negative in (|3|). For < k < 1, (|l|) refers to 
a ID periodic wavetrain of kinks. The kink sizes in the 
wavetrain are all equal to C^VT+T?, and the wavetrain 
period is: 



(14) 



where K(k)(and E(k) hereafter) is the Gomplete Jacobi 
Elliptic Integral. In the limit k — > 0, the kink wavetrain 
decays into harmonic oscillations whereas in the opposite 
limit i.e.K 1, it collapses turning to a single kink. The 
single kink describes a localized spatial density hole in 
the ground state of the repulsive BEG system. Remark 
that in this last limit the period L+ becomes infinite. 
It is interesting to calculate the number of particles in 
this kink wavetrain. We note from (|ll) that the wave- 
train function vanishes at a; = while being maximum, 
i.e. ip'^'^{x) = V'o^(k) at each period L+. In virtue of 
this periodicity we can readily shift the position of the 
wavetrain so that its density becomes maximum at its 
centre of mass. By this shift not only we avoid processes 
involving unsaturable particle densities, but also we pro- 
vide the wavetrain with a finite-energy rest frame. This 
leads to the following expression of the particle number: 
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Ki = 1 — K 



(15) 



This particle number in the kink wavetrain carries a total 
energy amounting to: 



1 + K- 



;K{k)]{p-pof'^ (16) 
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This is precisely the amount of energy needed to create 
the periodic wavetrain of kink condensate. 
Similarly, for < n < 1, ( |T3| ) will describe a ID periodic 
wavetrain of pulse(envelope) structures, the mean sizes 
of which are C,~\/2k^ — 1. The period of this envelope 
wavetrain is: 



L_ = ^2^2 - li^(fc)C+ (17) 
and its particle number and energy are respectively: 



E-{k) 



2h? 



[E{k) - nlKik)] 

(18) 
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-K(k) + 
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+6n^E{k)]{^io-^if^^ (19) 

As the first, this solution also reduces to harmonic os- 
cillations in the limit k — s- but transforms to a single 
envelope soliton as k — > 1. This single envelope soliton 
describes a localized spatial density hump forming in the 
ground state of the system in the attractive regime. 
We now turn to equation (^ describing the spatial dress- 
ing. We set ip'^ = ip'^f'^{x), where /'^(x) = f^^{x) and 

fonctions equation (|^) becomes 



^j^. In terms of the dimensionless 



Fiz) = 



(20) 



This is nothing but the eigenvalue equation ([^) with n = 
2. As functions of the physical parameters in the present 
model, we have: 
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P+(«,Pi) = (l-^i^)(l-|-«2), 



P„(k,p_l) = 4^2 + 1 
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(22) 



In terms of these parameters, the bound state solutions 
of (^ ) follow from the method presented above. In both 
cases they read by turns: 



/(3 4)(2^) = A3,4)(«:) sn^{z) - 
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P(3,4)(k, Pi) = 2(1 + ^2)^ 
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Tables | and |l^ list exact expressions of A^'' and p'"^}'^ ■ 
Parameters in these tables are calculated using the or- 
tho normalization relation for the Jacobi Elliptic func- 
tions l2^. We find: 
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flo = (1 + k'^)E{k) - kIK{k) 

ai = (1 + kI)E{k) - 2kIK{k) 
as = (2^2 - 1)E{k) + k\K{k) 
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(28) 

(29) 
(30) 



2, 6[5K{K)Ti{K{K)- E{k))± 
±{1 + k^)K{k)\, 5^ = 1-k^k\ (31) 



We end with a discussion on the asymptotic behaviours 
of the essential parameters calculated above when n = 1. 
For the repulsive regime, the kink wavetrain solution ( |Tl| ) 
will become: 



il)i'+{x) = V'^+(l)tanh 



VV2Cx^ 



l 9Po 
(32) 



and the associate bound states reduce to: 



fl^(x) ^ A,il)sech\^), pI^^O (33) 
/r(x) = A,il)sech{^-) tanh(-^). 



f^{x) = f^{x), pI,=pI,, 



-3^2 



2C 



(+)2 



(34) 
(35) 



As for the attractive case, the envelope-soliton wavetrain 
obtained in (13) becomes: 



r~{x) = ro-{i)sech (^^^ , ^r(i) = 

while the associate bound states reduce to: 



P^{x) = A,{l)sech\ — ), = 



20" 



(-)2 



/r(x)=Ai(l)secM — )tanh(— ), 
f^{x)=f^{x), p1i-p12 = 



(36) 

(37) 

(38) 
(39) 
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The meaning of these results is that for each regime there 
is at least one transverse fluctuation mode that couples 
to the single-soliton without energy exchange. This zero- 
energy mode then contributes in increasing the soliton 
stability, costing only a shape correction via uniform 
translation of the soliton centre of mass. The second 
bound state of the repulsive regime has negative energy 
and implies cost of energy to the soliton, which sees its 
creation energy lowered thus becoming unstable. On the 
contrary, the non-zero bound-state mode of the attrac- 
tive regime is positive and then contributes to the sta- 
bility of the bunch. In fact these two distinct behaviours 
can indicate which of the two regimes is properly de- 
scribed within the IDNS equation. From what precedes, 
the repulsive GP equation gives and unstable ID spatial 
soliton phase and therefore is unadequate for the study 
of the spatial localized density-hole excitations (dark soli- 
tons) in the ground state of the QID repulsive BEG 
system. It is instructive remarking that several other 
authors recently raised the same question |30 3^. The 
values of the particle numbers in this limit also can be 
helpful for checking the consistency of the theory. In the 
repulsive case the quantity ( |T5| ) gives iV+(l) = Podx, 
where (Ik = "sf^Ct is the actual width of the single kink 
resulting from As for the attractive case expres- 

sion (|19|) reduces to iV-(l) = 2poCx ■ The fact that the 
number of particles in the one-soliton condensate limit is 
the product of the condensate atomic density time the 
soliton size, agrees with the characteristic length scales 
defined above: rescaling the x axis in units of the coher- 
ence lengths gives the soliton densities (but no more 
the particle number in the soliton) po and 2po for the two 
respective regimes. 
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TABLE I. Amplitudes A'^q i 2 3 4 of the five bound states: 
= vT+"k? and = \/2k? — 1. The a, are defined in the 
text. 



regime 








^(3,4) 


repulsive 
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2 








4a2C^d+ 


2a(3,4) 


attractive 








TT 











Ci12 

TABLE IL Eigenvalues p\ 1234 of ^'^^ bound 
states. 













regime 


P±,o 




P±,2 


p1(3,4) 



repulsive ^(^+)2^^(+)2 ^(+)2^r+)2 ^^^2 [l ± 2c((+)^(«)] 
attractive ^(-f2^1-)2 ^r^r^^T^ [l T 2d(^(«)] 
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